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A covariant quark model, based both on the spectator formahsm and on vector meson dominance, 
and previously calibrated by the physical data, is here extended to the unphysical region of the lattice 
data by means of one single extra adjustable parameter - the constituent quark mass in the chiral 
limit. We calculated the Nucleon (A'^) and the yN — > A form factors in the universe of values for 
that parameter described by quenched lattice QCD. A qualitative description of the Nucleon and 
"/N — » A form factors lattice data is achieved for light pions. 



I. INTRODUCTION 

In principle QCD is the fundamental theory that de- 
scribes hadronic systems, but its non-perturbative char- 
acter makes a direct application to hadrons difficult, with 
the exception of the high energy and momentum trans- 
fer regime. Fortunately, at low effective field theo- 
ries based on effective hadronic degrees of freedom have 
been applied with success. Examples are chiral per- 
turbation theory [H, small scale expansion B, soliton 
models and constituent quark models [1, UlHl- Con- 
stituent quark models provided consistent descriptions 
encompassing the low-energy baryon spectrum, elastic 
and inelastic form factors, charge and magnetic radii, 
magnetic moments, axial and pseudoscalar form factors 

Although all those different frameworks may include 
the essential features of QCD, which are confinement and 
chiral symmetry, they are only partial simulations of the 
true underlying theory. 

Recently, significant progress has occurred in calcula- 
tions of QCD in the lattice, which evaluate the important 
QCD non-perturbative contributions at low directly 
from the underlying theory. Up to now the applications 
of lattice QCD are restricted still to large pion masses 
(tott > 350 MeV) corresponding to heavy quarks, and 
lattice spacings that are one order of magnitude smaller 
than the size of the nucleon (a ~ 0.05 fm). To extract 
information on the real world the results must be extrap- 
olated both to the continuum limit (a — *■ 0) and to the 
physical pion mass regime P, [l3|. However, except for 
the limit to the physical region, there is no simple way of 
interpreting the lattice QCD data. 

In this work we invert this procedure. We take here 
the challenge of describing the lattice data using a quark 
model. We start with the constituent quark model ob- 
tained from the covariant spectator theory. Our model 
was presented in Refs. 0, where it was adjusted to 
the experimental data for the four nucleon elastic form 
factors and the dominant form factor of the •jN — > A 
transition. In this work we extend this model to the 



unphysical region of the lattice data. Although the con- 
stituent quark models were originally thought for and 
applied to the physical limit, a constituent quark model 
that is simultaneously consistent with the experimental 
data and lattice QCD data is valuable, because it in- 
cludes indirectly the constraints of QCD, and therefore 
satisfies properties of the underlying theory. A similar 
procedure was considered in Refs. [la , [l^ using a larger 
number of parameters. 

Our model does not include explicit pion cloud ef- 
fects. However, as the electromagnetic interaction with 
the quarks is parametrized according to vector meson 
dominance (VMD), part of the pion cloud effects are in- 
directly taken into account. Therefore the model can be 
applied in the regions where the valence quark degrees of 
freedom are dominant, and it gives a good description of 
the nucleon and A systems [ij, [13 ■ 

For the nucleon, pion cloud effects are fundamental in 
the time-like region [I^l , although they are not so impor- 
tant in the region of small positive values [l3| , which 
is involved in the study of the electromagnetic transi- 
tion. For the — > A transition, as it happens also 
to other constituent quark models, without explicit pion 
cloud effects, the predictions for the magnetic moment, 
which dominates the transition, deviate necessarily from 
the data for low Q^, showing that an effective pion cloud 
within the VMD current is not sufficient. Because of 
the opening of the ttN channel, an explicit pion cloud 
contribution must be considered 0, [^. Therefore our 
results in the physical low region give only the con- 
stituent quark core contributions Q, usually labeled as 
"bare" contributions. Our results for these contributions 
are consistent with independent calculations of dynami- 
cal reaction models [l^ based on hadronic degrees of 
freedom. As for the E2 and C2 multipoles, in the large 
Nc limit ^] they represent second order corrections and 
are not considered here. 

As pion cloud effects in lattice QCD are expected to 
be small for > 0.40 GeV [2^], the description of the 
lattice QCD data appears as the ideal testing ground 
for our constituent quark model, where pion cloud ef- 
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fects are not explicitly included. This defines the main 
goal of this work. At present, lattice calculations are 
still restricted to high pion masses > 350 MeV, al- 
though technical improvements are increasingly allowing 
to reach nearer and nearer the physical region psi . [23 | . In 
this work we took the quenched lattice QCD data from 
Refs. [13, nil- In the unquenched calculations the effects 
of the sea quarks are explicitly considered. Although un- 
quenched lattice QCD data is already available, it is not 
expected that the quenched and unquenched data differ 
substantially in the region m^r > 0.40 GeV [22]. In ad- 
dition, there are some differences between unquenched 
data for the 77V A transition based on two differ- 
ent unquenched methods at similar pion masses (26| that 
still have to be clarified in the future. For these reasons 
we took the conservative option of using quenched data 
exclusively. 

Our new model presented here is based on three spe- 
cific features: i) the electromagnetic interaction with the 
constituent quarks is described within the impulse ap- 
proximation, and considering VMD; ii) the wavefunctions 
of the quark-diquark system are parametrized by simple 
monopole factors reduced to the Hulthcn form, with one 
or two effective range parameters that balance the de- 
tails of the short range and the long range behavior of 
the system; iii) the constituent quark magnetic anoma- 
lous moment scales with the inverse of its mass, which we 
write, following (27j as a function of the current quarks 
mass. 

A current based on VMD is suitable to describe the 
interaction of the photon with the constituent quark, for 
the quark-antiquark spin-1 vertex. Depending on the 
isospin, the intermediate meson pole corresponds to the p 
or the u! meson, at low Q^, or, in the large regime, to 
some other effective heavy meson pole Mh- The photon 
interaction is then described as proceeding through the 
production of an intermediate meson state which annihi- 
lates subsequently into a quark-antiquark pair. VMD is 
successful in the description of the electromagnetic inter- 
action with nucleons. 



II. EXTENSION OF THE QUARK MODEL TO 
THE LATTICE DATA REGION 

The nucleon S-state wavefunction includes the correct 
spin structure for the quark-diquark spin and 1 compo- 
nents, associated to isospin and 1 states, respectively. 
For the A, since total isospin is 3/2, the S-state wave- 
function reduces to the diquark spin 1 with isospin 1 
structure. As for the scalar wavefunction, it takes the 
phenomcnological form 



Nf 



ms{ai + XB){a2 + Xb)"" ' 



(1) 



where B = N,A, xb = ^^'^'"'m!^'^''''' . with M the 
baryon mass (for the nucleon or A), is the diquark 



mass and Nb is a normalization constant. [In Ref. [1, [31 
we considered ns = 1 for the nucleon and ub = 2 for the 
A]. The parameters ai can be interpreted as Yukawa 
mass or range coefficients that distinguish between two 
different regimes for the momentum range. 

The parametrization of the momentum dependence 
in terms of xs absorbs the dependence on the baryon 
masses M . The range parameters ai for the nucleon and 
the A were fixed by the nucleon and 77V A form fac- 
tor experimental data [1, [l3|. In the spectator quark 
model the transition amplitude between a initial (mo- 
mentum PJ) nucleon (N) and a final (momentum P+) 
baryon B = A can be written, in impulse approxima- 
tion [1, [1, [13 as 



= 3 



(2) 



where 4* b {P, k) is the complete baryon wavefunction (in- 
cluding spin, isospin and momentum dependence), A is 
the diquark polarization and jj the quark current oper- 
ator dependent of the hadron isospin /. The baryon po- 
larizations are suppressed for simplicity. The symbol J^. 
represents the invariant integral in the on-shell diquark 
moment = J 2E^{2Try-^ ' ^i^h Eg as the diquark energy. 
The factor 3 takes account for the isospin symmetriza- 
tion. 

The quark current [1, [^ [ij] takes the general form 



j7 



g/2+(Q')+2/2-(Q'K3 



2M 



N 



(3) 



The constituent quark form factors fi± {i = 1,2) are 
parametrized using a VMD structure, and are normalized 
according to /i±(0) = 1 and /2±(0) — k±, where k± 
are defined in terms of the quark anomalous magnetic 
moments k„ and according with = 2k„ — Kd and 
K- = |(2k„ + Kd). See Refs. 0, [31 for details. 

In this first calculation the wavefunctions are reduced 
to S-statcs for the quark-diquark system. Although it is 
well known that angular momentum components besides 
the S-states are essential for the description of the nu- 
cleon in Deep Inelastic Scattering, their effects are not 
so evident in the elastic form factors, particularly for low 
[l3|. As for the A, calculations of the valence quark 
contribution associated with higher angular momentum 
states (D-states) suggest a small effect ^5^ 2M- In fact, 
even with only S-wave components in the wavefunctions 
for both the nucleon and the A, the model generates the 
dominant contribution for the 7A'^ — > A transition, the 
dipolc magnetic moment form factor. The non-zero con- 
tributions for the other form factors appear only when D 
states are included, or when pion cloud effects are taken 
in consideration [i,[2^. At the end, we will estimate the 
effect in the — > A quadrupole transition form fac- 
tors, in the lattice regime, from adding the D-states to 
the present model. 
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As shown in fli*! the diquark mass scales out from the 
formulas obtained for the electromagnetic form factors. 
This allowed us to ignore any explicit dependence on the 
quark mass. This mass dependence was present in the 
nucleon, A and vector meson masses in an implicit way 
only. However, in order to grant a comparison of our 
results with the lattice data, we extend here the model 
to include a dependence of the quark anomalous moment 
K„ and Kd on the quark mass. This dependence was not 
explicitly considered in Refs. [1, [ij] because they dealt 
only with the physical data. 

Inspired by Ref . pTj , which combines chiral symmetry 
with conventional quark models, and applies an analytic 
continuation of the chiral expansion to the simple SU(6) 
model, we then use the smooth variation of the hadronic 
properties with the current quark masses above 60 MeV. 
Therefore, in the spirit of the constituent quark models, 
we assume here that Kq {q = u,d) scales with 1/Mq, 
where Mq is the constituent quark mass. Labeling the 
quark anomalous moment at the physical point by kJ^, 
we can write Kq, for an arbitrary constituent quark mass 



Mq, as 



Mq 



(4) 



where M^^'^ is the constituent quark mass at the physical 
point. 

To include an explicit dependence on the quark mass 
we consider the parametrization due to Cloet et. al. j27j 



Mq = M^ + cmq = M^ + cmf ^— (5) 

where mq (rn^^) is the (physical) current quark mass, 
is a new parameter corresponding to the constituent 
quark mass in the chiral limit (vriq = 0), and c is a coeffi- 
cient of the order of the unity. In the same equation 
stands for the pion mass in the model, a parameter to be 
fixed by the lattice data, and rri^^ for the physical mass. 
Following Ref. [23] we considered cmf'^ = 5.9 MeV. 

The parametrization in ^ is most sensitive to M^. 
Reference [27| fixes M^ = 0.42 GeV. Different descrip- 
tions using quark models and different lattice sizes lead 
to different values [15]. For this reason we use M-^ as 
the only free parameter allowed to vary in the calcula- 
tion presented here. It is needed to introduce an ex- 
plicit dependence of the nucleon magnetic moment on 
the pion mass, and to enable the connection of the con- 
stituent quark model to the lattice QCD calculations 

[Ulllilii,!!!,!!!. 

We consider three different cases. First we considered 



the case M^ 



~oo, corresponding to the limit where Kq 



has no dependence on m^. We tested also the original 
parametrization M^ = 0.42 GeV [13 • Finally, although 
we did not performed a systematic fit, we tested several 
other values of M^ . 



Furthermore, to extend the spectator model to the re- 
gion of the quenched lattice QCD data we still need to 
consider the nucleon and A masses determined by the 
lattice simulations [IJ] • This brings to the calculation an 
implicit dependence on the pion mass. As for the vecto- 
rial mesons included in the VMD, quark current picture, 
we use the parametrization (28| 



Co 



(6) 



where co = 0.776 GeV and ci = 0.427 GeV-^. The 
simple parametrization ^ describes well the quenched 
lattice QGD data and is consistent with finite volume 
corrections [Tsl . [30| . In this applications we consider this 
parametrization together with model II of the Refs. [1, 
Il4| . where the heavy vectorial meson mass is Mh = '2Mn, 
with Mn the nucleon mass of the lattice calculations. 



III. RESULTS 

In Fig. [1] we compare the predictions of our model 
for the nucleon form factor with the lattice data from 
Ref. [2^, corresponding to the three values of M^. We 
consider in particular the isovector nucleon form factor 
because the contributions of the disconnected diagrams 
vanishes in lattice calculations, if the flavor SU(2) sym- 
metry is assumed as in Ref. 



23| 



We considered the 
lowest pion masses and the smallest lattice spacings. In 
particular, we take the data corresponding to to^ — 504 
MeV, 649 MeV (a = 0.059 fm) and = 615 MeV 
(a — 0.078 fm). For larger lattice spacings a dependence 
on a is observed [23]. From the figure we conclude that 
M^ — 0.42 GeV and 0.80 GeV gives an excellent descrip- 
tion of the nucleon data. 

Still, when the pion mass increases there is a systematic 
deviation of our model from the lattice data of Ref. [1^ . 
This deviation may be a consequence of the fact that the 
wavefunction parametrization in terms of a low momen- 
tum scale (a^) (long range behavior) and high momen- 
tum scale (a^) was kept unchanged. As the pion masses 
vary, at least the long range parameter may vary, and 
may vary more rapidly for larger pion masses. As for 
the short range parameter (a^ ) it sets the scale below 
which the short range physics becomes important. As 
observed in applications of the finite-range regularization 
effective field theory, the scale associated to the momen- 
tum cut-off (and consequently short range effects) can 
be expressed by an universal regulator which describes 
simultaneously the larg e pion mass regime and the phys- 
ical regime [3,111 Ef. In our calculation is related 
to that universal regulator, and it is expected not to de- 
pend crucially on the pion mass value. 

For a finer analysis of our results, we compare our pre- 
dictions for the isovector nucleon magnetic moment in 
physical nucleon magnetons, with the lattice QCD data 
and the chiral result from Ref. [s^]- To convert G](^(0) 
given by the lattice data in units {^n nucleon 
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mass in lattice) to "physical" units ^^fpisr, we need to 

use the transformation Gli{0) ^^G'Xf(O). The re- 
sults are presented in the Fig. [2] for M^^ = 0.42 GeV and 
= 0.80 GeV, as function of m^. The agreement of our 
results with the chiral expression from Ref. shows the 
consistency of our calculations with chiral calculations. 

Finally, in Fig.[3]we compare the quenched lattice data 
for the ■fN A transition from Ref. [l^l with the spec- 
tator quark model corresponding to = 0.42 GeV, 
0.80 GeV and — +oo. All cases shown correspond 
to a lattice spacing of a = 0.092 fm. The contribution of 
the quark core extracted from ^2^ at the physical point 



is also included. To be consistent, we considered the 
parametrization of Eq. ^ for nip, although the origi- 
nal Ref. [24j gives slightly different results. As for Af^r 
and Ma , we use the values derived directly from the lat- 
tice data . For the larger pion masses we observe 
an almost perfect agreement between the predictions of 
Af^ = -t-oo and the data, although Af^ = 0.80 GeV is 
also close. For = 0.411 GeV we have also a good 
agreement, except for a slight deviation from the lattice 
data for 0.7 GeV^ < < 2 GeV^. This deviation may 
result from the effect of the pion cloud for light pions, 
predicted to be important for < 0.40 GeV [22|. As 
for the physical pion mass case, our model is coherent 
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FIG. 2: Nucleon isovector magnetic moment as function of 
nin, in physical nucleon magnetons. Theoretical line from 
Ref. [s^l- Our results are labeled "VMD" and shown for two 
values of the parameter. 



with the constituent quark core data labeled 'Bare Form 
Factor', which is extracted indirectly from experiment 
0,123] ■ In this case all lines coincide. 

As mentioned already, the D-states in the A wave- 
function induce contributions for the subleading electric 
and Coulomb quadrupole transition form factors and 
G^. The exact contributions depend on the specific 
parametrization, in particular on the admixture coeffi- 
cients for the two D-waves. In Ref. [l^ it was shown 
how a percentage of D-states of « 0.9% can provide an 
excellent description of the quadrupole lattice data, with- 
out significantly changing the description of the domi- 
nant magnetic dipole form factor. That application cor- 
responds to the = +00 limit and estimates the ef- 
fect of the valence quark (or bare) contributions as only 
w 20% of the total quadrupole for both G^ and G^, 
at the physical point (the remaining contribution being 
the pion cloud). On the other hand, in that particular 
parametrization, and in the region of study < 1.5 
GeV^, the effect of the constituent quark mass depen- 
dence described by Eq. ^ corresponds to a correction of 
less than 20% for both form factors, a correction smaller 
than the typical lattice errorbands (« 30%). 



IV. CONCLUSIONS 

Chiral perturbation methods for lattice QCD extrapo- 
lations are useful for small (like < 1.5 GeV^), but 
are not adequate for the high region. Constituent 
quark models can supply an alternative guidance for lat- 
tice QCD extrapolations. In this work we consider a co- 
variant constituent quark model of the nucleon and the A 
based on the spectator formalism for the quark-diquark 



system (covariance is an important issue for the descrip- 
tion of the high momentum transfer processes). 

The model contains no explicit pion cloud, besides the 
effects included in the p term of the VDM picture for the 
electromagnetic current, and is therefore reduced to the 
bare quark hadron structure. Since pion cloud effects are 
expected to be negligible for large values of the pion mass, 
the comparison of the model to the quenched data for 
pion masses larger than 450 MeV is justified a priori. To 
accomplish that comparison the initial covariant quark 
model was extended to the lattice data region, m,r < 700 
MeV. In light of the work in Ref. [l^l , this extension was 
done by introducing the constituent quark mass in the 
chiral limit parameter, M^, and the nucleon p, A masses 
used in the lattice calculations. 

The main conclusion is that the covariant constituent 
quark model which was previously calibrated by means 
of a quantitative description of the nucleon, — > A 
and A form factor data in the physical region, as shown 
in previous works [1, [13, after a simple extension 
involving one parameter only, describes also quantita- 
tively well the lattice results for the nucleon isovector 
form factor and the — > A "bare" magnetic form fac- 
tor. With = 0.42 GeV, consistently with the range 
Af^ = 0.3 — 0.45 GeV suggested by several constituent 
quark models [1, 0, H, we obtain a very good de- 
scription of the lattice nucleon form factors data, but 
underestimate the "/N —^ A magnetic moment form fac- 
tor data by less than 9% at low for to^ < 500 MeV. 
Note that the original parametrization — 0.42 GeV 
[2?! was a result of a phenomenological fit and was not 
derived from first principles. 

An optimal description of the lattice data for both nu- 
cleon and 7-/V — > A transition form factors is achieved, 
once the scale of the constituent quark mass in the 
chiral limit is fixed as = 0.80 GeV. [The values 
= 0.42 and 0.80 GeV are better for the nucleon data; 
= 0.80 GeV and +00 are better for the ■^N ^ A 
data]. The parameter M^^ — 0.80 GeV is relatively large 
when compared with alternative constituent quark mod- 
els, which may indicate that some of the effective quark- 
antiquark configurations contained in our model through 
the VMD mechanism in the electromagnetic current, for 
instance, do not correspond to the lattice calculations 
in the quenched approximation. Still, the possibility of 
adjusting other parameters is very promising. 

A refinement of the description can be pursued, by in- 
creasing the number of adjustable parameters, or simply 
by introducing an extra dependence on the pion mass. 
For instance, in the future we certainly plan to check 
the assumption that the ai parameter which controls 
the long range regime does not depend on the pion mass 
value. In addition, we may consider as well an explicit 
dependence of the heavier meson mass pole Mh, that reg- 
ulates the short range effect in VMD mechanism, on the 
pion mass. Alternatively, a pion mass dependence of the 
a2 parameter can be introduced. Furthermore, we want 
to study the quality of the description in the high pion 
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masses region, not probed yet here. 
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